Ultrawide phononic band gap for combined in-plane and out-of-plane waves 
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We consider two-dimensional phononic crystals formed from silicon and voids, and present opti- 
mized unit cell designs for (1) out-of-plane, (2) in-plane and (3) combined out-of-plane and in-plane 
elastic wave propagation. To feasibly search through an excessively large design space (~10 40 pos- 
sible realizations) we develop a specialized genetic algorithm and utilize it in conjunction with the 
reduced Bloch mode expansion method for fast band structure calculations. Focusing on high- 
symmetry plain-strain square lattices, we report unit cell designs exhibiting record values of nor- 
malized band-gap size for all three categories. For the combined polarizations case, we reveal a 
design with a normalized band-gap size exceeding 60%. 



Phononic crystals (PnCs) are periodic materials that 
exhibit distinct frequency characteristics such as the pos- 
sibility of formation of band gaps. Within a band gap, 
wave propagation is effectively prohibited. This inher- 
ent dynamical phenomenon can be utilized in a broad 
range of technologies at different length scales. Appli- 
cations of PnCs include elastic/acoustic waveguiding [lj 
and focusing vibration minimization [3J, sound colli- 
mation [I] , frequency sensing [5] 15] , acoustic cloaking [7J , 
acoustic rectification |Hj , opto- mechanical waves coupling 
in photonic devices [9], thermal conductivity lowering in 
semiconductors [TUHTB"] . among others [Ti] . 

In general, it is most advantagous to have the fre- 
quency range of a band gap maximized while pulling its 
midpoint as low as possible in order to keep the unit cell 
size to a minimum. Selecting the topological distribu- 
tion of the material phases inside the unit cell provides a 
a powerful means towards reaching this target, and this 
has been the focus of numerous research studies not only 
on PnCs but also photonic crystals (PtCs). 

The exploration for optimal unit cell designs was ini- 
tiated by Cox and Dobson in 1999 [15] (in the context 
of PtCs). The articles by Burger et al. [H] and Jensen 
and Sigmund [17] provide a review of subsequent studies 
concerned with band-gap widening in PtCs. In the area 
of PnCs, the problem has been treated in a variety of 
settings and using several techniques. For example, unit 
cells have been optimized in one-dimension |18U19| and in 
two-dimensions (2D) [20H2S], using gradient-based [2"TT - 
[23] as well as non-gradient-based [24] [25] techniques. In- 
terest in band-gap size maximation has also been treated 
outside the scope of the unit cell dispersion problem 
[2T1 126] . In all these optimization studies the focus has 
been primarily on PnCs based on an infinite thickness 
model and a material composition consisting of two or 
more solid (or solid and fluid) phases with the exception 
of a few investigations that considered thin-plate single- 
phase models [25] [53J . Recognizing the practical signif- 
icance of solid-and-air PnCs with relatively large cross- 
sectional thickness, some studies considered the configu- 
raton of a 2D solid matrix with periodic cylindrical voids 
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- modeled under 2D plain-strain conditons |27] or as a 
three-dimensional continuum with free surface bound- 
ary conditions [28] - and investigated the dependence 
of band-gap size upon the void radius. For combined 
out-of-plane and in-plane waves in 2D infinite-thickness 
PnCs formed from silicon and a square lattice of circular 
voids, it has been shown that the band-gap size normal- 
ized with respect to the mid-gap frequency cannot exceed 
40% [27] . In this letter we utlize a specialized optimiza- 
tion algorithm in pursuit of the best unit cell solid-void 
distribution for the 2D plain-strain problem considering 
high-symmetry square lattices. We cover the cases of 
(1) out-of-plane, (2) in-plane and (3) combined out-of- 
plane and in-plane elastic wave propagation. Our search 
methodology is also applicable to the parallel problem 
of 2D PtCs optimization, where transverse-electric and 
transverse-magnetic waves may be considered separately 
[29] or in combination [30] . 

The governing continuum equation of motion for a het- 
erogeneous medium is 

V.C: i(Vu+(u) T )-pu, (1) 

where C is the elasticity tensor, p is the density, u is the 
displacement vector, x = {x,y, z} is the position vector, 
V is the gradient operator, and (.) T is the transpose op- 
eration. We assume the wave propagation to be confined 
to the x-y plane only, that is, du/dz = 0. As such, we 
have two independent sets of equations, one for out-of- 
plane motion and the other for in-plane motion. To ob- 
tain the band structure for a given PnC unit cell design 
we assume a Bloch solution to the governing equations 
in the form u(x,k;t) = u(x, k)e l ( k ' x ~ wt ) where u is the 
Bloch displacement vector, k is the wave vector, u> is the 
frequency, t is the time, and i = y— 1. Due to lattice 
symmetry the analysis is restricted to the first Brillouin 
zone. We consider square lattices and furthermore im- 
pose C± v symmetry at the unit cell level. Subsequently 
design representation is needed in only a portion of the 
unit cell and the band structure calculation is limited to 
the corresponding irreducible Brillouin zone (IBZ). Fur- 
themore, we model only the solid portion of the unit cell. 
The void portion is not modeled since we permit only 
contiguous distribution of solid material. In this manner 
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the PnCs considered exhibit geometric periodicity (with 
free in-plane surfaces) and not material periodicity. In 
practice the voids will be either in vacuum or filled with 
air. Our model presents an adequate representation of 
both cases because the elastic waves propagating in the 
solid will have the dominating effect [28]. In fact, this also 
suggests that the results we show are practically indepen- 
dent of the choice of the solid material. We numerically 
solve the emerging eigenvalue problem using the finite 
element (FE) method utilizing 4-node bilinear quadrilat- 
eral elements. In addition, the reduced Bloch mode ex- 
pansion (RBME) method [3T] is applied to the FE model 
to substantially speed up the band structure calculations 
throughout the optimization process. In the RBME im- 
plementation we use a two-point expansion. The final 
reported results however are based on full (non-reduced) 
calculations. 

Unit Cell Optimization We represent a square unit cell 
Y by nxn pixels forming a binary matrix G. This ma- 
trix is then reduced in size following the underlying unit 
cell symmetry. Each of the pixels can be assigned to 
either a no-material (void) or a material (silicon), i.e., 
cjij G {0, 1}. Throughout all the intermediate steps of the 
optimization process, we treat the void pixels as a highly 
compliant medium and this enables us to conveniently 
manipulate the unit cell designs. Once the optimization 
is complete, we assess the final designs by modeling only 
the silicon portion of the unit cell as described above. 

The objective function is formulated in terms of the 
size of a particular band-gap width normalized with re- 
spect to its midpoint frequency: 

_ maxjmin^iujf^ik^g)) - mox^^ffe.g)), 0) 

m - (m»;! 1 « 1 (^j)) + ™<ik 2 (^s)))/2 

(2) 

where min^^uf {kj , g)) and maXjjLi(u)f (kj , g)) denote 
the minimum and maximum, respectively, of the i th 
frequency uji over the entire discrete wave vector set, 
kj,j = 1, ... ,rifc, tracing the border of the IBZ. The band 

gap exists only when the minimum of the (i + l) th branch 
is greater than the maximum of the i th branch; otherwise 
no band gap exists. 

We employ a genetic algorithm (GA) to maximize f(g). 
A GA is a nature-inspired optimization technique that 
mimics biological evolution. It generally starts with a 
pool of candidate solutions (i.e., designs) according to a 
certain objective (fitness) function, then applies a group 
of operators, namely, selection, crossover and mutation, 
in order to evolve to more fit designs (i.e., with higher 
objective function values). Compared to gradient-based 
methods, GAs are less likely to get trapped into local 
minima especially for problems with vast search spaces, 
hyper dimensions and large number of variables |32j . as 
is the case in our unit-cell optimization problem. 

In our GA the initial population of unit-cell designs 
is set up to be random to avoid any initial bias that 
might negatively affect the search. Since it is unlikely to 
have a band gap at the onset, the area between the two 



dispersion branches of interest is subsequently used as an 
indicator of the fitness of the unit cell design: 

Fitness = H a + <t>\f {jg) + foHi + <t> 3 f c {g), (3) 

where subscript (.)i. 0jC denotes the wave type (i.e., either 
out-of-plane, in-plane or combined), ^1-3 are constants 
equal to 10 4 , 10 8 , 10 15 respectively, introduced to set pri- 
orities during the evolution process, and Hj is a step 
function defined as: 

H _ f if f(g) > (band gap exists) ,^ 
j \ Area if f (g) = (no band gap) ^ ' 

In Eq. Q the Area represents a measure of the "area" 
in frequency-wavenumber space between the two consec- 
utive dispersion branches of interest: 

Area = £ [(^ +1 (k Jl9 )) - (w 2 (%, <?))] (5) 
i=i 

The only condition enforced in the initial population is 
that adjacent pixels of the same material type appear in 
pairs in each row. Throughout the evolution, tournament 
selection and single-point crossover are the two types of 
operations applied on any given pair of "parent" unit-cell 
designs. Following the unit-cell symmetry constraint, the 
"offspring" mutates according to a specific probability 
using the rule: Select random pixel x; if 53r=-l 9x+r > 1, 
set each of the three pixels to one, otherwise, to zero. 
The GA terminates when no further improvement in the 
objective function value is noted for a prescribed number 
of generations. At the end of the search, the final unit- 
cell topology passes through a simple one-point flip local 
search for fine tuning and smoothening. 

Lead-Follow Algorithm The combined out-of-plane and 
in-plane optimization problem poses a challenge in set- 
ting up the objective function because it is based on two 
sets of independent equations. Here we adopt a unique 
strategy, which we refer to as a lead-follow algorithm, 
whereby the search for a combined band gap is tackled in 
a two-stage fashion during the evolution process. The al- 
gorithm starts with a set of random designs and searches 
for a band gap for out-of-plane waves (the leader) be- 
tween two prescribed branches guided by Eq. ^ for 
indication of design quality. Once the GA opens a gap, 
it shifts its focus to the in-plane waves (the follower), 
but now the branch numbers encompasing the band gap 
are not prescribed - they are determined by the same 
frequency range that spans the band gap of the leader 
wave type. The objective function for the follower wave- 
type is set to be also the Area as given by Eq. This 
lead-follow process continues until a combined band gap 
is found, at which point the objective function effectively 
switches to being the actual value of the normalized com- 
bined band gap. This process is automated through a 
generalized fitness function as given in Eq. (|2|. We note 
that in principle the identity of the leader and the fol- 
lower can be reversed. 
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Results In applying the specialized GA we considered 
the following properties for isotropic silicon (A and /j, de- 
note Lame's coefficients): p s = 2330 Kg/m 3 , A s = 85.502 
GPa, /i s = 72.835 GPa, and we used a resolution of 
ri = 32. At this resolution, the total number of possi- 
ble material distributions within the unit cell domain is 
8.7xl0 40 . This highlights the tremendously large search 
space that the GA needs to navigate through. At the end 
of each complete GA run, we doubled the resolution of 
the emerged topology to become 64x64 pixels, and then 
smoothened the topology (while keeping it pixelated) by 
following a few simple rules. For the third case (i.e., com- 
bined waves) a splines-based solid material distribution 
has been subsequently generated to represent a manufac- 
turable design. Figure 1 presents the unit cell topologies 
and band structures of the optimized unit cells for the 
three cases, and Table 1 lists the objective values ob- 
tained. In our results we identify a band gap number by 
the number of the optical branch that borders it from the 
top. 



TABLE I. Normalized band-gap (BG) size for the optimized 
unit cells 



Wave Type Out-of-Plane In-Plane Combined 
BG Number 1st 2nd 2nd Lowest 
Representation Pixels Pixels Pixels Pixels Splines 
Norm. BG Size 1.2270 1.1132 0.7696 0.6259 0.6021 



The optimized unit cell topologies for out-of-plane 
waves (Figs. la,b) show contiguous solid media ap- 
proaching the limiting case of isolated square or circu- 
lar inclusions. This limiting case represents the optimal 
conditions for sonic crystals which admits only pressure 
waves |33j . For the presented case, the thin connections 
shown are needed to support the propagation of the shear 
elastic waves. The optimized topology for the in-plane 
waves problem on the other hand show a mostly solid 
material with delicately shaped voids. This is consistant 
with the understanding in the literature that solid mate- 
rial with isolated voids represents the optimal conditions 
for band-gap opening for in-plane waves |33] . We note 
that no band gap appeared below the first optical branch 
due to the difficulty in preventing this branch from cross- 
ing through the acoustic branches. The optimal design 
for the combined case appears to be a blend (although 
non-intuitive in shape) among the out-of-plane and in- 
plane design traits. Upon appropriate size scaling to the 
frequency range of interest, all designs are amenable to 
fabrication by splines-based smoothening with minimal 
loss in objective value (as demonstrated in Fig. le). 
When compared to corresponding solid-void PnC unit- 
cell configurations reported in the literature, each of the 
designs shown in Fig. 1 represent a record value of nor- 
malized band-gap size in its category. 
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FIG. 1. Optimized unit cell design and band structure for 
out-of plane waves: (a) 1st band gap, (b) 2nd band gap, in- 
plane waves: (c) 2nd band gap, combined out-of-plane and in- 
plane waves, lowest band gap: (d) pixels and (e) splines. The 
minimum feature size and radius of curvature are identified in 
(e). /* denotes the objective function modified to represent 
the normalized band-gap size for the combined waves case. 
All band gaps are shaded in grey. 
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